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Chapter  1 
Introduction 


Accurate  knowledge  of  the  sea  surface  directional  spectrum  is  important  for  oceano¬ 
graphic  and  remote  sensing  applications.  For  example,  the  sea  surface  directional 
spectrum  is  used  in  the  composite  surface  model  to  predict  sea  surface  backscattering, 
and  various  sea  surface  directional  spectral  models  may  yield  very  different  backscat¬ 
tering  results.  It  is  the  purpose  of  this  report  to  evaluate  sea  surface  backscattering  in 
the  microwave  region  for  three  surface  directional  spectra  using  the  composite  surface 
method. 

The  composite  surface  model  has  been  proven  to  be  a  powerful  and  efficient  tool 
for  obtaining  useful  quantitative  sea  surface  backscattering  results  [l]-[5] .  Under  the 
composite  surface  (or  two-scale)  model,  the  ocean  surface  spectrum  is  separated  into 
long  and  short  wave  regions  (relative  to  the  electromagnetic  wavelength)  through 
choice  of  a  separation  wavenumber,  typically  chosen  as  2  to  3  electromagnetic  wave¬ 
lengths.  Near  normal  incidence  radar  returns  are  described  in  this  method  according 
to  geometrical  optics,  while  oblique  incidence  radar  cross  sections  are  obtained  by 
averaging  small  perturbation  method  predictions  over  the  slope  distribution  of  the 
long  wave  spectral  region.  The  method  has  limitations  as  near  grazing  incidence  ge¬ 
ometries  are  approached,  but  should  provide  accurate  predictions  for  grazing  angles 
approximately  ten  degrees  or  larger  [6]. 

The  three  sea  surface  wavenumber  spectral  models  investigated  in  this  report  are 
a  power  law  spectrum,  the  Durden- Vesecky  spectrum  [7]  and  the  Apel  spectrum  [8]. 
The  power  law  spectrum  describes  a  sea  surface  which  does  not  vary  with  the  wind 
velocity  over  the  sea  surface.  Although  it  is  not  a  realistic  description  of  the  sea 
surface,  its  simplicity  is  utilized  to  demonstrate  the  viability  of  the  composite  surface 
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method  to  predict  the  sea  surface  backscattering.  This  is  done  by  comparing  results 
obtained  from  the  composite  surface  method  with  those  obtained  from  the  small  slope 
approximation  for  rough  surface  scattering  [9].  The  small  slope  approximation  is  a 
completely  different  method  from  the  composite  surface  method  for  calculating  sea 
surface  scattering.  The  former  does  not  require  a  separation  of  the  spectrum  into 
long  and  short  wave  regions,  but  it  is  not  as  efficient  as  the  latter  in  terms  of  the 
amount  of  computations  required.  Because  they  are  two  completely  different  methods, 
agreement  between  results  obtained  by  them  should  indicate  the  applicabilities  of  both 
methods. 

Unlike  the  power  law  spectrum,  the  Durden- Vesecky  and  Apel  spectra  do  vary  with 
the  wind  velocity  over  the  sea  surface.  The  differences  between  the  latter  two  spectra 
lie  primarily  in  the  short  wave  portion  of  the  spectrum  and  the  azimuthal  variations 
of  the  spectrum.  One  glaring  difference  between  them  is  observed  in  the  fact  that  the 
Durden- Vesecky  model  places  all  azimuthal  variations  in  short  gravity-capillary  and 
capillary  waves  only,  while  the  Apel  model  places  azimuthal  variations  in  both  long 
and  short  waves.  Thus  the  Durden- Vesecky  and  Apel  spectra  are  good  candidates  for 
studying  sea  surface  directional  spectral  model  influence  on  sea  surface  backscattering. 
Because  they  describe  more  realistic  sea  surfaces,  backscattering  results  obtained  for 
these  two  spectra  can  be  compared  with  empirical  scattering  models  available  in  the 
literature.  Expressions  for  an  additional  surface  spectral  model  recently  developed  by 
the  Naval  Research  Laboratory  [10]-[11]  are  also  included,  but  scattering  calculations 
with  this  model  have  not  yet  been  completed. 

Chapter  2  of  this  report  considers  the  expressions  needed  for  application  of  the 
composite  surface  method  and  reviews  the  surface  spectral  models  considered.  Sea 
surface  backscattering  results  obtained  from  the  composite  surface  method  with  the 
power  law  spectrum  are  then  compared  with  those  obtained  from  the  small  slope  ap¬ 
proximation  in  Chapter  3.  Comparisons  of  backscattering  results  at  1  GHz,  5  GHz, 
10  GHz,  14  GHz  and  35  GHz  for  the  Durden- Vesecky  and  Apel  spectra  with  wind 
speeds  ranging  from  4  m/sec  to  20  m/sec  at  a  height  of  19.5  m  above  the  sea  surface 
are  given  in  Chapter  4.  Backscattering  results  at  5  GHz  and  14  GHz  are  also  com¬ 
pared  with  three  empirical  scattering  models  (CMOD2-I3,  SASS-II,  and  NSCAT-I) 
in  Chapter  4. 
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Chapter  2 

The  Composite  Surface  Method 
And  Sea  Surface  Spectral  Models 


This  chapter  briefly  reviews  the  formulation  of  the  composite  surface  method  for 
calculating  the  backscattered  radar  cross  sections  from  a  sea  surface,  and  also  provides 
expressions  for  the  sea  surface  spectral  models  considered  in  the  study. 

2.1  Composite  Surface  Method 

As  mentioned  in  the  previous  chapter,  the  composite  surface  method  separates  the 
sea  surface  wavenumber  spectrum  into  long  and  short  wave  regions  (relative  to  the 
electromagnetic  wavelength).  In  this  report,  the  separation  wavenumber  is  chosen  to 
be  one  half  of  the  incident  electromagnetic  wavenumber,  which  has  been  shown  to 
yield  favorable  scattering  results  compared  to  those  obtained  by  a  rigorous  numerical 
method  [12].  Near  normal  incidence,  the  long  wave  spectral  region  yields  specular 
returns  according  to  geometrical  optics,  while  oblique  incidence  radar  returns  are 
obtained  by  averaging  small  perturbation  method  predictions  (Bragg  scattering  from 
the  short  wave  spectral  region)  over  the  slope  distribution  of  the  long  wave  spectral 
region. 

Let  W(kx,ky )  denote  the  sea  surface  wavenumber  spectrum.  For  simplicity,  the 
a;-axis  and  y-axis  are  in  the  up-wind  and  cross-wind  directions,  respectively,  with  the 
z-axis  pointing  upward  from  the  sea  surface.  The  slope  variances  of  the  sea  surface  in 
the  up- wind  and  cross-wind  directions  due  to  its  long  wave  spectral  region  are  given, 
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respectively,  by 


s2  =  /  dkp  /  d(j>  k3p  cos2  <p  W(kp  cos  <p,  kp  sin  (p)  (2.1) 

Jo  Jo 

r\k  r  27T 

Sy=  dkp  /  dtp  k3  sin2  <p  W ( kp  cos  </>,  kp  sin  (p)  (2.2) 

Jo  Jo 

where  k  is  the  wavenumber  of  the  incident  electromagnetic  wave,  kp  cos  <p  =  kx  and 
A;p  sin  <p  =  ky. 

The  specular  radar  return  due  to  the  long  wave  spectral  region  of  the  sea  surface 
spectrum,  assuming  a  Gaussian  slope  distribution,  is  given  by 


a  = 


irp 


—  |  tan2  6{ 


2  cos4  di  sxsy 


cos2  <f>i  ,  sin2  <f>j 
=  ?  '  ,2 


(2.3) 


where  a  is  the  normalized  radar  cross  section,  ( Oi ,<pi)  are  the  incident  polar  and 
azimuthal  angles  with  9i  measured  from  the  z-axis,  and  |T|2  is  the  power  reflection 
coefficient  from  a  calm  sea  surface  at  normal  incidence, 

in2  =  Tt4.  2  (2-4) 


1  + 


with  e  being  the  relative  dielectric  constant  of  the  sea.  A  higher  order  correction  to 
T,  which  is  small,  has  been  omitted  in  (2.4).  The  expression  in  (2.3)  is  valid  for  both 
horizontal  and  vertical  polarizations.  It  is  seen  from  (2.3)  that  the  specular  radar 
return  drops  quickly  as  the  incident  angle  9i  increases.  Thus  (2.3)  is  significant  only 
for  near  normal  incidence. 

The  expressions  for  the  Bragg  scattering  radar  cross  sections  from  a  tilted  sea 
surface  are  considered  next.  Define  a  primed  coordinate  system  such  that  the  x'- 
and  y'- axes  are  in  and  perpendicular  to  the  plane  of  incidence,  respectively;  i.e.  the 
x'y'  coordinate  system  is  obtained  by  rotating  the  xy  coordinate  system  by  the  angle 
4>i  about  the  z  axis.  The  tilted  sea  surface  is  defined  by  two  angles,  ip  and  S.  The 
angle  ip  is  obtained  by  rotating  the  untilted  sea  surface  normal  (z-  or  z'-axis)  about 
the  y'-axis  resulting  in  a  double  primed  (x"y"z")  coordinate  system,  and  the  angle 
5  is  obtained  by  rotating  the  new  surface  normal  (z"-axis)  about  the  x"-axis.  The 
incident  angle  with  respect  to  the  tilted  surface  becomes  0[,  which  is  related  to  9{  by 


cos#'  =  cos(9i  +  ip)  cos 8. 


(2.5) 
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Now,  according  to  the  small  perturbation  method,  the  Bragg  scattering  radar  cross 
sections  are  given  by  [13] 


—  167T&4  COS4#' 
° vv(0'i)  =  167t/c4  cos4  #' 


a  cos  5  \ 2 
sin  O'i  ) 


a  cos  6  \ 2 
sin#'  ) 


W(k'x,k'y )  (2.6) 
W(k'x,k'y)  (2.7) 


kL 


e  —  1 

(cos  9[  +  yJT—  sin2  #')2 


_ €-_l _ 

(ecos#(  +  y/T—  sin2  #')2 
2k  (a  cos  </>i  —  7  sin  $  sin  fa 


[e(l  +  sin2  0') 

) 


-  sin2  #'] 


k'y  =  2k  (a  sin  fa  +  7  sin  5  cos  fa) 
a  =  sin(0j  + 1/)) 


7  =  cos(0j  + 1/0 


(2.8) 

(2.9) 

(2.10) 

(2.11) 

(2.12) 

(2.13) 


The  expressions  in  (2.6)  and  (2.7)  are,  respectively,  copol  normalized  radar  cross 
sections  for  the  horizontal  and  vertical  polarizations.  Note  that  W(k'x,ky )  in  (2.6) 
and  (2.7)  should  be  set  to  zero  if  (k'x  +  k'y)  <  \k  as  in  the  cases  of  near  normal 
incidence,  which  means  that  (2.6)  and  (2.7)  are  significant  only  for  oblique  incidence 
and  are  the  contributions  of  the  short  wave  spectral  region  of  the  sea  surface  spectrum. 
These  radar  returns,  (2.6)  and  (2.7),  should  be  averaged  over  the  slope  distribution  of 
the  long  wave  spectral  region.  Assuming  a  Gaussian  slope  distribution,  the  averaged 
radar  cross  sections  for  fa  =  0  or  7r/2  are 


\  0Vv(0i) 


}  ■  LxdztLjzy{°t 


) 


27r  SXSy 


(  £2L 

V  si 


(2.14) 


where  (zx  =  tan  ip,zy  =  tan<5)  for  fa  =  0  and  ( zx  =  tan  5,zy  —  t&mp)  for  fa  =  7t/2. 
The  above  integration  is  performed  numerically  using  Gauss-Legendre  quadrature. 
Because  of  the  Gaussian  distribution,  the  integration  ranges  in  (2.14)  can  be  truncated 
to  \zx\  <  4 sx  and  \zy\  <  4 sy. 

For  the  cross  polarized  radar  return  calculation,  a  strict  application  of  the  com¬ 
posite  surface  method  would  require  adding  first  and  second  order  small  perturbation 
results  for  a  tilted  sea  surface  averaged  over  the  long  wave  slope  distribution.  These 
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calculations  would  require  a  prohibitive  amount  of  computation.  Instead,  a  simple 
second  order  small  perturbation  solution  for  an  untilted  sea  surface  as  obtained  by 
Valenzuela  [14]  will  be  used  as  follows: 


where 


&Hv(0i)  —  87t/c4COS 2  6i 


/oo  pc 

dfox  I 

•OO  J  — c 


2(*?  -  k2fkuikz 


(ku,  +  ka)(k\k,i  +  k2kui) 


dk „ 


( kxiky  kyj  kX)(  kx{  kX  kyj  ky  ) 


(kli  +  KXklh  +  Vku) 

IV ( k%  kxi)  ky  kyi)  fV (k>x  fcy  4”  ^ yi ) 


k\  —  ek2 
kXi  =  A:sin0jCos</>i, 


kyi  =  k  sin  9i  sin  4>i 


kzi  —  k  cos  Qi 


hzi  =  ik\  -  k2xi  -  k2yif12 

kz  =  (k2  -  k2x  -  kl)1/2 
hz  =  (kj  -  k2x  -  k2y)1/2 

With  the  substitution  of  kx  —  kpcos(j)  and  ky  =  kpsin4 >,  (2.15)  becomes 

(e  —  l)2y/e  —  sin2  6>, 


(2.15) 


(2.16) 

(2.17) 

(2.18) 

(2.19) 

(2.20) 


&Hv{Qi)  =  87rA:4cos4^i 


(cos  9i  +  \ft  —  sin2  9t)  (e  cos  9t  +  \fe  —  sin2  9t) 


f00  JL  f 2T  Sln2  “  &) 

■y.  "Jo  ^  Iet  +  *1,I2 

(^x  ^xi>  ky  (^x  ”1”  ^ yi ) 


(2.21) 


Note  the  integrals  above  are  over  the  entire  surface  spectrum  so  that  both  long  and 
short  wave  contributions  are  included;  a  validation  of  the  above  procedure  will  be 
provided  in  Chapter  3  through  comparison  with  the  small  slope  approximation.  Again 
the  required  numerical  integration  is  performed  using  Gauss-Legendre  quadrature. 


2.2  Sea  Surface  Spectral  Models 

The  formulation  of  the  composite  surface  method  given  above  will  be  used  to  evaluate 
radar  returns  from  a  rough  sea  surface  for  three  sea  surface  wavenumber  spectral 
models:  a  power  law  model,  the  Durden- Vesecky  model  and  the  Apel  model.  The 
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power  law  spectral  model  is  described  by 

W ( kp  cos  <f>,  kp  sin  <j>) 


a0/kj 


(2.22) 


where  oo  =  This  simple  spectral  model  does  not  vary  with  the  wind  velocity  or 
azimuthal  angle  and  is  primarily  used  to  validate  the  composite  surface  method  by 
comparing  scattering  results  with  those  obtained  from  the  small  slope  approximation. 

The  Durden- Vesecky  spectrum  was  derived  in  [7]  to  match  sea  surface  backscat- 
tered  data  under  a  composite  surface  model,  and  thus  is  not  derived  from  direct 
hydrodynamic  measurements.  The  spectrum  is  described  by 


W (kp  cos  </>,  kp  sin  (f)  = 


-[1  +  C  (1  —  e~ss'kp)  cos  20] 


.0.74 

\“19.5 kP  J 


2\  alogio(*p/2) 


The  various  parameters  in  (2.23)  are  defined  as  follows: 

a\  —  0.004, 
ss  =  1.5  x  10-4, 


for  kp  <2 
for  kp  >  2 


(2.23) 


a  =  0.225, 

g  =  9.81  +  7.25  x  10'5A:J, 

=  2(1  -  R)  1 

1  +  R  1  -  gT 
^  0.003  +  0.00192ui2.5 

0.00316wi2.5  ’ 

/o"  dkt  k‘Akt) 

—  kpW(kp  cos  <j>,  kp  sin  <^) 

P*  ~  1  +  C(1  —  e_S5'fc?)  cos2(/) ' 

Also  note  that  n*  is  the  surface  friction  velocity  in  m/s  and  the  wind  velocity  uh  at 
a  height  h  (in  meters)  above  the  sea  surface  is  related  to  u*  by 


0.4  L  6.84  x  10 -5/u*  +  4.28  x  lO"3^  -  4.43  x  10-4J 
The  Apel  spectrum  was  derived  in  [8]  to  match  measured  buoy  long  wave  spectra 
and  measured  wave  tank  short  wave  spectra,  and  thus  has  a  more  direct  connection 
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to  hydrodynamic  data.  However,  other  researchers  have  criticized  slope  moments  ob¬ 
tained  from  this  spectrum,  as  well  as  the  large  amplitudes  predicted  near  the  gravity¬ 
capillary  wave  phase  speed  minimum  [15].  The  model  is  described  by 


kp 


L  ^e-(V^-V'M2/0.32fcp1 


10[— 4.95+3.45(1- 


0.8 kp  sech 


kp  -  400 


1  +  (Aip/100)2 

,  for  0  <  4>  <  7t/2 


io/4-7)] 

(2.24) 


where  a 2  =  0.00195  and  kp  —  •  Note  that  (2.24)  is  valid  for  0  <  (f>  <  7t/2. 

The  spectrum  in  the  other  quadrants  is  generated  by  reflecting  (2.24)  about  the  x- 
and  y-axes.  The  factor  A  in  (2.24)  is  introduced  to  keep  the  surface  height  variance 
constant.  A  is  given  by 


7Ta/ 0.14  +  5(kp/kp)1-9  /erf 


where  erf  is  the  error  function. 

A  final  model  developed  from  hydrodynamic  measurements  by  the  Naval  Research 
Laboratory  [10]- [11]  is  also  considered  and  will  be  included  in  future  comparisons 
of  scattering  measurements.  The  NRL  model  provides  a  one-dimensional  surface 
spectrum  which  can  be  extended  to  two  dimensions  if  a  form  for  surface  azimuthal 
variations  is  assumed.  If  a  uniform  cos2  (f>  variation  is  assumed  for  simplicity,  the 
resulting  two  dimensional  spectrum  is 

4  <  *p  <  16 

2.01x  10~2 
1.19  % 

1.97x10-*  l-3  1 

1.19  g+rk 2 

2,Ti1Q4  V 6  900  <  kP  <  1600 

where  k  is  in  rads/m ,  g  is  the  acceleration  of  gravity,  and  r  is  the  ratio  of  water 
surface  tension  to  its  density.  The  above  form  is  based  on  tests  which  insure  that  an 
integration  of  the  two  dimensional  spectrum  over  ky  approximates  the  original  one¬ 
dimensional  NRL  spectrum.  Inclusion  of  recently  proposed  azimuthal  distributions 
for  this  spectrum  [16]-[17]  are  currently  under  consideration. 


16  <  kp  <  100 

100  <  kp  <  900  ^  ’ 
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2.2.1  Comparison  of  Spectral  Models 

To  compare  spectra  predicted  by  the  Durden- Vesecky,  Apel,  and  NRL  models,  ap¬ 
proximate  the  directional  spectrum  by 

W ( kp  cos  4 b,  kp  sin  <j>)  «  Wo(kp)  +  W?,{kp)  cos  2 <f>.  (2.26) 

Figures  2. 1-2.3  show  the  comparisons  of  the  corresponding  zeroth  and  second  har¬ 
monic  curvature  spectra,  i.e. ,  co(kp)  =  kpWo(kp)  and  C2(p)  =  kpW2(kp),  at  three  wind 
speeds.  While  the  three  models  show  agreement  in  the  zeroth  harmonic  to  within  ap¬ 
proximately  an  order  of  magnitude,  there  are  significant  differences.  The  NRL  model 
is  observed  to  predict  the  largest  amplitudes  in  the  10  <  k  <  100  range  for  the  higher 
wind  speeds,  and  to  have  a  peak  for  k  >  100  which  is  similar  to  that  of  the  Apel  spec¬ 
trum.  Both  the  NRL  and  Apel  models  typically  predict  larger  spectral  amplitudes 
than  the  Durden- Vesecky  model,  except  for  very  small  or  very  large  wavenumbers. 
The  second  harmonic  term  is  responsible  for  azimuthal  variation  in  the  spectrum. 
One  can  see  that  the  Durden- Vesecky  model  places  all  azimuthal  variations  in  the 
short  wave  region,  while  the  Apel  model  has  azimuthal  variations  in  both  long  and 
short  wave  regions.  Second  harmonics  for  the  NRL  model  are  not  plotted  due  to  the 
simple  cos2  <f>  assumption  used. 

The  directional  spectrum  at  the  Bragg  scattering  wavenumber  is  important  as  it 
influences  oblique  incidence  radar  cross  sections.  The  Bragg  scattering  wavenumber 
is  2fcsin0j,  which  depends  on  the  incident  frequency  and  incident  angle.  Figure  2.4 
shows  the  comparisons  of  W(2A;sin7r/4, 0)  for  the  Durden- Vesecky,  Apel,  and  NRL 
spectral  models  at  five  frequencies.  These  comparisons  show  the  NRL  model  to  pre¬ 
dict  significantly  larger  Bragg  wavenumber  amplitudes  for  the  frequencies  considered, 
except  at  the  highest  frequency  where  Apel  spectrum  results  are  similar.  The  Apel 
model  amplitudes  are  larger  than  those  of  the  Durden- Vesecky  model,  with  the  great¬ 
est  degree  of  similarity  of  5  GHz.  These  results  suggest  that  copol  radar  returns  at 
oblique  incidence  will  be  significantly  larger  with  the  NRL  spectral  model  than  with 
the  Apel  and  Durden- Vesecky  models,  and  somewhat  larger  for  the  Apel  model  than 
the  Durden- Vesecky  model.  Detailed  comparisons  of  scattering  cross  sections  for  the 
Apel  and  Durden- Vesecky  spectra  are  presented  in  Chapter  4. 
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Zeroth  harmonic  at  4  m/s  Second  harmonic  at  4  m/s 


Wavenumber  (rads/m)  Wavenumber  (rads/m) 

Figure  2.1:  Comparisons  of  curvature  spectra  at  U19.5  =  4  m/sec. 


Zeroth  harmonic  at  1 2  m/s  Second  harmonic  at  1 2  m/s 


Wavenumber  (rads/m)  Wavenumber  (rads/m) 

Figure  2.2:  Comparisons  of  curvature  spectra  at  «19.5  =  12  m/sec. 
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Spectral  amp  (m4)  Spectral  amp  (m4)  Spectral  amp  (m4)  <§’  Curvature  Spectrum  Amplitude 


Zeroth  harmonic  at  20  m/s  Second  harmonic  at  20  m/s 


ire  2.3:  Comparisons  of  curvature  spectra  at  U19.5  =  20  m/sec 


io-10> — i - : - ; - 1 

5  10  15  20 


Windspeed  (m/s) 


Windspeed  (m/s) 


Figure  2.4:  Comparisons  of  fF(2A;sin0;,O)  for  0*  =  7t/4 


Chapter  3 

Radar  Cross  Sections  for  the 
Power  Law  Spectrum 


As  a  test  of  the  applicability  of  the  composite  surface  model,  backscattering  results 
are  compared  in  this  chapter  with  those  obtained  from  the  small  slope  approxima¬ 
tion  [9].  The  small  slope  approximation  (SSA)  is  a  more  rigorous  theory  of  rough 
surface  scattering  which  can  include  both  large  and  small  scale  roughness  without 
choice  of  an  arbitrary  separation  wavenumber  in  the  spectrum.  The  SSA  is  based 
on  a  perturbation  expansion  in  surface  “quasi-slope”,  with  the  zeroth  order  term 
producing  an  expression  for  average  cross  sections  similar  to  that  obtained  from  the 
Kirchhoff  approximation  with  a  modified  polarization  dependence  outside  of  the  sur¬ 
face  integration.  Calculation  of  the  first  order  correction  to  zeroth  order  results 
requires  an  additional  integration  over  the  surface  Fourier  transform,  thus  increasing 
computational  requirements.  Because  SSA  expressions  are  more  difficult  to  evaluate 
than  those  of  the  composite  model,  use  of  the  composite  model  is  more  efficient  if 
it  can  be  shown  that  the  two  models  yield  similar  results.  The  relationship  between 
the  first  order  SSA  and  the  composite  model  has  been  considered  in  [18],  where  it  is 
shown  that  the  two  should  yield  similar  predictions  if  the  first  order  SSA  correction 
is  included. 

To  simplify  SSA  model  calculations,  a  power  law  spectrum  is  used  in  the  compar¬ 
ison.  Since  the  power  law  spectrum  does  not  vary  with  wind  speed  or  azimuth  angle, 
the  number  of  parameters  in  the  comparison  is  substantially  reduced.  Note  that 
the  simple  power  law  wavenumber  spectrum  given  in  equation  (2.22)  grows  without 
bound  as  kp  approaches  zero.  Thus  it  is  necessary  to  impose  a  low  cutoff  wavenum¬ 
ber  such  that  the  spectrum  exists  only  for  kp  >  kcntos.  Variations  in  backscattering 
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cross  sections  with  frequency  can  also  be  neglected  in  this  comparison  if  the  cutoff 
wavenumbers  chosen  are  described  relative  to  the  electromagnetic  wavenumber  k. 

To  avoid  difficulties  in  first  order  SSA  model  calculations,  a  Monte  Carlo  eval¬ 
uation  of  SSA  cross  sections  with  100  realizations  was  performed  following  the  de¬ 
terministic  surface  method  of  [19].  Since  the  deterministic  surfaces  generated  are  of 
finite  size,  the  maximum  surface  size  and  thus  minimum  low  cutoff  wavenumber  are 
limited  by  computational  capacity.  The  results  shown  used  512  by  512  point  surfaces 
of  size  128  by  128  electromagnetic  wavelengths.  Two  comparisons  using  low  cutoff 
wavenumbers  k/ 128  and  k/ 64  are  performed  to  examine  the  effect  of  increasing  sur¬ 
face  low  frequency  content.  For  surface  rms  height  h,  the  resulting  kh  products  are 
5.7  and  2.9  respectively  so  that  surfaces  can  be  considered  moderately  rough  on  an 
electromagnetic  scale.  Surface  permittivities  are  chosen  to  model  sea  water  at  14 
GHz. 

Backscattered  radar  cross  sections  for  low  cutoff  wavenumber  k/ 128  (2.25  rad/m 
at  14  GHz)  calculated  using  the  composite  surface  method  are  compared  with  those 
obtained  by  the  small  slope  approximation  in  Figure  3.1.  The  agreement  between 
the  two  is  very  good  for  all  observation  angles  in  co-polarized  results,  indicating  the 
applicability  of  the  composite  model.  Predictions  of  cross  polarized  radar  returns 
obtained  using  equation  (2.21)  also  match  SSA  results  well,  except  for  6i  <  10°  where 
SSA  results  are  over  predicted.  Similar  conclusions  are  also  observed  when  the  low 
cutoff  wavenumber  is  set  at  k/ 64  (4.5  rad/m  at  14  GHz)  as  shown  in  Figure  3.2. 
Comparison  of  Figures  3.1  and  3.2  shows  the  expected  effect  of  increasing  surface  low 
frequency  content:  cross  sections  decrease  near  normal  incidence  while  increasing  at 
larger  observation  angles.  Overall,  these  results  indicate  the  applicability  of  both  the 
composite  and  SSA  methods  for  calculating  the  radar  cross  sections  from  a  rough 
sea  surface,  and  motivate  use  of  the  composite  model  for  larger  scale  studies  of  sea 
surface  backscattering. 
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freq=1 4  GHz,  power  law  spectrum 


theta  angle  (deg) 


Figure  3.1:  Comparisons  of  radar  returns  for  the  power  law  spectrum  with  low  cutoff 
wavenumber  kj  128. 


freq=14GHz,  power  law  spectrum 


Figure  3.2:  Comparisons  of  radar  returns  for  the  power  law  spectrum  with  low  cutoff 
wavenumber  k/ 64. 
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Chapter  4 

Radar  Cross  Sections  for  the 
Durden- Vesecky  and  Apel  Spectra 


The  Durden- Vesecky  and  Apel  spectral  models  vary  with  the  wind  velocity  over  the 
sea  surface  and  thus  describe  more  realistic  sea  surfaces.  To  evaluate  the  effect  dif¬ 
fering  models  for  the  sea  surface  spectrum  can  have  on  sea  surface  scattering,  radar 
cross  sections  for  the  two  spectral  models  are  calculated  using  the  composite  surface 
method  at  five  frequencies  (1  GHz,  5  GHz,  10  GHz,  14  GHz  and  35  GHz)  and  for 
three  wind  speeds  (C/19.5  of  4  m/sec,  12  m/sec  and  20  m/sec).  Results  are  illustrated 
for  upwind  cross  sections  and  the  up-cross  wind  difference,  and  also  compared  with 
predictions  of  empirical  backscattering  models. 

4.1  Upwind  cross  sections 

Figure  4.1-4.15  compare  predictions  of  backscattering  under  the  composite  model  for 
the  Durden- Vesecky  and  Apel  spectra  versus  observation  angle.  Results  at  L-band 
(Figures  4. 1-4. 3  for  wind  speeds  4,  12,  and  20  m/s,  respectively)  show  significant 
differences  in  predictions  from  the  two  spectra,  with  the  Apel  spectrum  generally 
predicting  smaller  and  larger  cross  sections  for  normal  and  oblique  observation,  re¬ 
spectively.  Larger  oblique  observation  cross  sections  for  the  Apel  model  are  consistent 
with  the  comparison  of  Bragg  wavenumber  amplitudes  presented  in  Chapter  2.  Sim¬ 
ilar  results  are  obtained  at  C-band  (Figures  4. 4-4. 6),  X-band  (Figures  4. 7-4.9),  Ku- 
band  (Figures  4.10-4.12),  and  Ka  band  (Figures  4.13-4.15).  The  two  spectra  generally 
are  in  closest  agreement  at  the  lowest  wind  speed  considered,  with  the  exception  of 
L-band,  and  produce  the  smallest  differences  for  other  wind  speeds  at  C-band.  Again, 
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these  results  are  consistent  with  the  comparison  of  Bragg  wavenumber  amplitudes  in 
Chapter  2. 

Variations  with  frequency  show  somewhat  smaller  cross  sections  near  normal  inci¬ 
dence  as  frequency  increases  and  somewhat  larger  cross  sections  at  oblique  angles,  as 
should  also  be  expected  for  surfaces  which  appear  rougher  for  shorter  electromagnetic 
wavelengths.  Cross  polarized  cross  sections  are  observed  to  increase  significantly  with 
wind  speed  and  slightly  with  frequency.  Note  that  the  cross  polarized  results  should 
be  inaccurate  for  0*  <  10°  as  discussed  in  Chapter  3. 

4.1.1  Normal  incidence  cross  sections 

Of  some  special  interest  to  altimeter  applications  are  radar  cross  sections  at  normal 
incidence.  As  indicated  in  equation  (2.3),  the  rms  slopes  sx  and  sy  are  important 
for  determining  near  normal  incidence  radar  returns.  Figure  4.16  compares  sx  and 
sjsxsy  associated  with  the  Durden- Vesecky  and  Apel  spectral  models  for  three  wind 
speeds  at  five  frequencies.  One  can  see  that  sx  and  are  always  larger  for 

the  Apel  model  than  for  the  Durden- Vesecky  model.  Thus,  from  equation  (2.3)  and 
Figure  4.16  one  can  reach  the  conclusion  that  near  normal  incidence,  the  radar  cross 
section  for  the  Apel  model  is  lower  than  for  the  Durden- Vesecky  model,  and  as  the 
incident  angle  increases  away  form  normal  incidence,  “tilting”  effects  will  be  larger 
for  the  Apel  model.  Normal  incidence  returns  from  Figures  4.1  to  4.15  are  plotted  in 
Figure  4.17  for  both  spectral  models  at  the  5  frequencies  and  3  wind  speeds  considered. 
Radar  cross  sections  at  normal  incidence  are  indeed  higher  for  the  Durden- Vesecky 
spectrum  than  for  the  Apel  spectrum,  and  results  for  both  spectral  models  tend  to 
decrease  as  the  wind  speed  picks  up  from  4  m/sec  and/or  the  frequency  increases 
from  1  GHz  to  35  GHz.  Note  the  frequency  dependence  for  these  near  specular 
results  is  obtained  from  the  choice  of  k/2  as  the  long  wave  region  cutoff  wavenumber, 
since  slope  variances  are  obtained  only  from  the  long  wave  region.  Comparisons  with 
empirical  altimeter  results  would  be  required  to  assess  these  predictions  in  detail. 
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Figure  4.1:  Comparisons  of  radar  cross  sections  in  the  up- wind  direction  at  1  GHz 
with  wind  speed=4  m/sec. 
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Figure  4.2:  Comparisons  of  radar  cross  sections  in  the  up-wind  direction  at  1  GHz 
with  wind  speed=12  m/sec. 
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Figure  4.3:  Comparisons  of  radar  cross  sections  in  the  up-wind  direction  at  1  GHz 
with  wind  speed=20  m/sec. 
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freq=5  GHz,  wind  speed=4  m/sec 


theta  angle  (deg) 

Figure  4.4:  Comparisons  of  radar  cross  sections  in  the  up-wind  direction  at  5  GHz 
with  wind  speed=4  m/sec. 
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Figure  4.5:  Comparisons  of  radar  cross  sections  in  the  up-wind  direction  at  5  GHz 
with  wind  speed=12  m/sec. 
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Figure  4.6:  Comparisons  of  radar  cross  sections  in  the  up-wind  direction  at  5  GHz 
with  wind  speed=20  m/sec. 
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Figure  4.7:  Comparisons  of  radar  cross  sections  in  the  up-wind  direction  at  10  GHz 
with  wind  speed=4  m/sec. 
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Figure  4.8:  Comparisons  of  radar  cross  sections  in  the  up-wind  direction  at  10  GHz 
with  wind  speed=12  m/sec. 


freq=10  GHz,  wind  speed=20  m/sec 


theta  angle  (deg) 

Figure  4.9:  Comparisons  of  radar  cross  sections  in  the  up-wind  direction  at  10  GHz 
with  wind  speed=20  m/sec. 
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freq=14  GHz,  wind  speed=4  m/sec 


theta  angle  (deg) 

Figure  4.10:  Comparisons  of  radar  cross  sections  in  the  up-wind  direction  at  14  GHz 
with  wind  speed=4  m/sec. 


freq=14  GHz,  wind  speed=12  m/sec 


theta  angle  (deg) 

Figure  4.11:  Comparisons  of  radar  cross  sections  in  the  up-wind  direction  at  14  GHz 
with  wind  speed=12  m/sec. 


freq=14  GHz,  wind  speed=20  m/sec 


theta  angle  (deg) 

Figure  4.12:  Comparisons  of  radar  cross  sections  in  the  up-wind  direction  at  14  GHz 
with  wind  speed=20  m/sec. 


freq=35  GHz,  wind  speed=4  m/sec 


theta  angle  (deg) 

Figure  4.13:  Comparisons  of  radar  cross  sections  in  the  up-wind  direction  at  35  GHz 
with  wind  speed=4  m/sec. 


freq=35  GHz,  wind  speed=12  m/sec 
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theta  angle  (deg) 

Figure  4.14:  Comparisons  of  radar  cross  sections  in  the  up- wind  direction  at  35  GHz 
with  wind  speed=12  m/sec. 


freq=35  GHz,  wind  speed=20  m/sec 


0  15  30  45  60  75 

theta  angle  (deg) 

Figure  4.15:  Comparisons  of  radar  cross  sections  in  the  up-wind  direction  at  35  GHz 
with  wind  speed=20  m/sec. 
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Figure  4.16:  Comparisons  of  rms  slopes  for  the  Durden- Vesecky  and  Apel  spectra. 
The  top  figure  is  for  sx  and  the  bottom  figure  for  yj s (Note  that  the  curves  for 
the  Durden- Vesecky  spectrum  are  simply  straight  line  connections  of  the  calculated 
data  points  at  five  frequencies.) 


Figure  4.17:  Radar  cross  sections  at  normal  incidence  for  the  Durden- Vesecky  and 
Apel  spectra.  (Note  that  the  curves  for  the  Durden- Vesecky  spectrum  are  simply 
straight  line  connections  of  the  calculated  data  points  at  five  frequencies.) 
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(a) 


Figure  4.18:  Wind  exponents  versus  (Bragg)  wavenumber  obtained  directly  from 
spectral  models  and  obtained  from  Upwind  RCS  data  (a)  VV  (b)  HH 

4.1.2  Wind  speed  dependence 

The  wind  dependence  of  radar  cross  sections  is  also  of  interest,  and  is  typically  mod¬ 
eled  as  cr  =  aU\ 9  5  (a  not  in  dB)  where  a  and  b  are  constants  which  vary  with  frequency 
and  observation  angle.  Figure  4.18  plots  the  exponent  b  versus  Bragg  wavenumber  (a 
function  of  frequency  and  observation  angle)  derived  from  fits  to  upwind  RCS  data 
with  30°  <  9i  <  75°  at  wind  speeds  4,  8,  12,  16,  20,  and  24  m/sec  and  frequencies  1, 
5,  10,  14,  and  35  GHz.  Results  from  backscatter  data  under  the  two  spectral  models 
(symbols)  are  compared  to  those  obtained  directly  from  the  spectra  (curves).  Expo¬ 
nents  derived  from  VV  RCS  data  are  found  to  be  in  general  agreement  with  those 
directly  from  the  spectra,  indicating  the  small  influence  of  tilting  effects  on  VV  cross 
sections.  HH  results  show  larger  differences  since  tilting  effects  are  more  significant, 
particularly  at  the  higher  frequency  end  of  the  Bragg  region  for  a  given  frequency 
corresponding  to  the  larger  observation  angles.  The  two  spectral  models  predict  gen¬ 
erally  increasing  exponents  versus  wavenumber,  although  the  Durden  Vesecky  trend 
is  more  linear  while  the  Apel  model  shows  a  higher  order  dependency. 
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4.2  Variations  with  azimuth  angle 


Figures  4.1  to  4.15  plot  radar  cross  sections  when  the  incident  angle  fa  is  in  the  up 
wind  direction,  i.e.,  fa  =  0.  Usually  the  difference  in  the  radar  cross  sections  for  the 
cross  wind  direction  (fa  =  n/2)  and  for  the  up  wind  direction  is  small.  Thus,  instead 
of  showing  the  radar  cross  sections  for  fa  =  7t/2,  which  would  be  very  similar  to  those 
shown  in  Figures  4.1  to  4.15,  it  is  more  informative  to  show  the  difference  in  the  radar 
cross  sections  for  the  two  fa  angles.  To  do  this,  write  sea  surface  radar  cross  sections 
(not  in  decibels)  as 

a  =  A  +  B  cos  fa  +  C  cos  2fa ,  (4.1) 

where  A,  B  and  C  are  functions  of  the  frequency,  wind  speed  and  incident  angle 
6i .  The  first  harmonic  term  B  describes  an  up/down  wind  asymmetry  in  radar  cross 
sections,  and  can  be  obtained  only  with  a  non-Gaussian  random  process  model  of  the 
sea  surface.  For  simplicity,  first  harmonic  terms  are  neglected  in  this  report.  The 
second  harmonic  term  C  is  responsible  for  the  radar  cross  section  variation  in  the  up 
wind  and  cross  wind  directions.  Examples  of  2C  at  5  GHz  and  14  GHz  are  shown 
in  Figures  4.19  to  4.22.  The  second  harmonic  term  has  larger  values  for  0*  <  30° 
where  cross  sections  in  general  are  larger,  and  shows  variations  with  polarization 
similar  to  those  obtained  from  cross  section  zeroth  azimuthal  harmonics  A.  Cross 
section  second  harmonics  show  improved  agreement  between  the  two  spectral  models 
for  oblique  observation  angles,  but  near  normal  incidence  results  are  much  larger  for 
the  Apel  model  due  to  its  inclusion  of  azimuthal  variations  in  the  long  wave  portion 
of  the  spectrum. 

4.3  Comparison  with  empirical  scattering  models 

It  is  interesting  to  see  how  the  radar  cross  section  predictions  of  the  Durden- Vesecky 
and  Apel  spectral  models  compare  with  some  empirical  sea  surface  scattering  mod¬ 
els.  Three  representative  empirical  models  in  the  literature  are  the  CMOD2-I3  [20], 
SASS-II  [21]  and  NSCAT-I  [22]  models.  The  CMOD2-I3  model  applies  for  vertically 
polarized  radar  cross  sections  at  5  GHz;  comparisons  of  upwind  cross  sections  with 
the  composite  model  are  shown  in  Figures  4.23  to  4.25  for  wind  speeds  4,  12,  and 
20  m/sec,  respectively.  Note  the  CMOD2-I3  empirical  model  strictly  applies  for  inci- 
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freq=5  GHz,  wind  speed=4  m/sec 


Figure  4.19:  Differences  between  the  normalized  radar  cross  sections  in  the  upwind 
and  cross  wind  directions  at  5  GHz  with  wind  speed=4  m/sec. 


freq=5  GHz,  wind  speed=12  m/sec 


Figure  4.20:  Differences  between  the  normalized  radar  cross  sections  in  the  upwind 
and  cross  wind  directions  at  5  GHz  with  wind  speed=12  m/sec. 
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freq=1 4  GHz,  wind  speed=4  m/sec 


Figure  4.21:  Differences  between  the  normalized  radar  cross  sections  in  the  upwind 
and  cross  wind  directions  at  14  GHz  with  wind  speed=4  m/sec. 


freq=14  GHz,  wind  speed=12  m/sec 


Figure  4.22:  Differences  between  the  normalized  radar  cross  sections  in  the  upwind 
and  cross  wind  directions  at  14  GHz  with  wind  speed=12  m/sec. 
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dence  angles  0*  in  the  range  18°  <0{  <  58°  only  even  though  the  empirical  curves  are 
plotted  outside  this  range.  Results  show  excellent  agreement  between  the  empirical 
model  and  cross  sections  from  the  Apel  spectrum,  while  Durden- Vesecky  spectrum 
results  underpredict  the  CMOD2-I3  model.  Second  azimuthal  harmonic  results  are 
compared  in  Figures  4.26  to  4.27  for  wind  speeds  4  and  12  m/sec,  respectively,  and 
produce  similar  conclusions.  The  success  of  the  Apel  spectrum  in  matching  CMOD2- 
13  results  suggests  that  the  C-band  Bragg  portion  of  the  spectrum  is  being  modeled 
reasonably  well  by  the  Apel  spectrum,  although  complete  conclusions  in  this  regard 
are  difficult  to  obtain  due  to  the  interplay  of  Bragg  scatter  and  “tilting”  effects. 

The  SASS-II  and  NSCAT-I  empirical  models  provide  predictions  of  both  horizon¬ 
tal  and  vertical  co-pol  cross  sections  at  14  GHz  for  observation  angles  0  <  0*  <  60° 
and  15°  <  0j  <  65°,  respectively.  These  models  provide  similar  but  slightly  differing 
predictions  since  they  were  derived  from  observations  with  different  sensors.  Compar¬ 
isons  with  the  SASS-II  model  for  up- wind  cross  sections  at  4  and  12  m/s  are  provided 
in  Figures  4.28  to  4.29  and  for  up/cross  wind  differences  in  Figures  4.30  to  4.31. 
The  comparison  shows  empirical  predictions  generally  to  fall  between  those  from  the 
Durden  Vesecky  and  Apel  spectral  models,  although  vertically  polarized  cross  sec¬ 
tions  are  somewhat  closer  to  those  from  the  Durden  Vesecky  spectrum.  The  large 
up/cross  wind  differences  obtained  from  the  Apel  spectrum  for  small  incidence  angles 
are  not  predicted  by  the  SASS-II  model,  suggesting  that  some  of  the  Apel  spectrum 
azimuthal  description  in  the  long  wave  region  may  be  inaccurate.  Comparisons  with 
the  NSCAT-1  model  in  Figures  4.32  to  4.35  yield  similar  results. 
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freq=5  GHz,  wind  speed=4  m/sec 


Figure  4.23:  Vertical  polarization  radar  cross  sections  in  the  up  wind  direction  at 
5  GHz  with  wind  speed=4  m/sec. 


freq=5  GHz,  wind  speed=12  m/sec 


Figure  4.24:  Vertical  polarization  radar  cross  sections  in  the  up  wind  direction  at 
5  GHz  with  wind  speed=12  m/sec. 


freq=5  GHz,  wind  speed=20  m/sec 


Figure  4.25:  Vertical  polarization  radar  cross  sections  in  the  up  wind  direction  at 
5  GHz  with  wind  speed=20  m/sec. 
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freq=5  GHz,  wind  speed=4  m/sec 


Figure  4.26:  Differences  between  the  vertical  polarization  radar  cross  sections  in  the 
up- wind  and  cross-wind  directions  at  5  GHz  with  wind  speed=4  m/sec. 


freq=5  GHz,  wind  speed=12  m/sec 


Figure  4.27:  Differences  between  the  vertical  polarization  radar  cross  sections  in  the 
up-wind  and  cross-wind  directions  at  5  GHz  with  wind  speed=12  m/sec. 
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Figure  4.28:  Radar  cross  sections  in  the  up-wind  direction  at  14  GHz  with  wind 
speed=4  m/sec. 
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normalized  radar  cross  section  (dB) 


theta  angle  (deg) 


Figure  4.29:  Radar  cross  sections  in  the  up-wind  direction  at  14  GHz  with  wind 
speed=12  m/sec. 


second  harmonic  term  2C  (dB) 


freq=14  GHz,  wind  speed=4  m/sec 


Figure  4.30:  Differences  between  the  radar  cross  sections  in  the  up-wind  and  cross- 
wind  directions  at  14  GHz  with  wind  speed=4  m/sec. 
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freq=14  GHz,  wind  speed=4  m/sec 
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freq=14  GHz,  wind  speed=12  m/sec 
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Figure  4.35:  Differences  between  the  radar  cross  sections  in  the  up-wind  and  cross- 
wind  directions  at  14  GHz  with  wind  speed=12  m/sec. 
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Chapter  5 
Conclusions 


The  composite  surface  method  has  been  employed  to  evaluate  sea  surface  backscat- 
tering  of  three  ocean  surface  spectral  models.  Radar  cross  sections  for  the  power 
law  spectrum  were  found  to  agree  well  with  those  calculated  by  the  small  slope 
approximation,  validating  the  composite  surface  model  for  prediction  of  scattering 
from  ocean-like  surfaces.  The  Durden- Vesecky  and  Apel  spectral  models  were  then 
applied  to  describe  more  realistic  sea  surfaces  which  vary  in  azimuth  and  with  wind 
speed.  These  two  models  were  found  however  to  differ  in  their  description  of  the  short 
gravity-capillary  portion  of  the  spectrum  and  in  surface  azimuthal  dependency.  Radar 
cross  sections  calculated  from  the  two  models  generally  showed  large  differences,  in¬ 
dicating  the  importance  of  the  surface  spectral  model  used  for  sea  surface  scattering 
predictions.  Comparisons  with  empirical  scattering  models  showed  the  Apel  model 
to  provide  better  agreement  for  C-band  results,  while  empirical  predictions  generally 
were  between  those  from  the  two  spectral  models  at  Ku  band.  Further  comparisons 
with  a  larger  set  of  empirical  results  can  help  to  clarify  the  successes  and  failures  of 
these  spectral  models,  and  will  be  considered  in  future  efforts.  Future  studies  will 
also  include  additional  directional  spectrum  models  to  evaluate  their  performance, 
and  will  investigate  sea  surface  thermal  emission  effects  as  well. 
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